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ABSTRACT 

We present a detailed analysis of 108 helium-line (DB) white dwarfs based on 
model atmosphere fits to high signal-to-noise optical spectroscopy. We derive a 
mean mass of 0.67 Mq for our sample, with a dispersion of only 0.09 Mq. White 
dwarfs also showing hydrogen lines, the DBA stars, comprise 44% of our sample, 
and their mass distribution appears similar to that of DB stars. As in our previous 
investigation, we find no evidence for the existence of low-mass (M < 0.5 Mq) 
DB white dwarfs. We derive a luminosity function based on a subset of DB white 
dwarfs identified in the Palomar-Green survey. We show that 20% of all white 
dwarfs in the temperature range of interest are DB stars, although the fraction 
drops to half this value above Teg ~ 20, 000 K. We also show that the persistence 
of DB stars with no hydrogen features at low temperatures is difficult to reconcile 
with a scenario involving accretion from the interstellar medium, often invoked 
to account for the observed hydrogen abundances in DBA stars. We present 
evidence for the existence of two different evolutionary channels that produce 
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DB white dwarfs: the standard model where DA stars are transformed into DB 
stars through the convective dilution of a thin hydrogen layer, and a second 
channel where DB stars retain a helium-atmosphere throughout their evolution. 
We finally demonstrate that the instability strip of pulsating V777 Her white 
dwarfs contains no nonvariables, if the hydrogen content of these stars is properly 
accounted for. 

Subject headings: stars: abundances — stars: evolution — stars: fundamental 
parameters — stars: luminosity function, mass function — stars: oscillations — 
white dwarfs 



Introduction 



DB white dwarfs are characterized by the extreme purity of their helium-dominated 
atmospheres with, occasionally, only minute traces of hydrogen (DBA stars) and heavy 
elements (DBQ or DBZ stars) seen in their optical spectra. It has been 15 years since 
we last revie wed the work carried ou t in Montreal on the optical spectra of DB and DBA 
white dwarfs (iBeaucharnp et al.lll996l). Parts of this work have been published over the years 



(|Beauchamp et al. 



1999 



Hunter et al. 



1995 



1997 



2001 



Beauchamp||l998uBeauchamp fc Wesemaellll998l : lBeauchamp et al 



Wesemael et al.ll200ll ). and full and definitive results from this long- 



term project are presented here. 

Progres s on several fronts relevant to the study of DB stars has been significant since 
our review (IBeauchamp et al.l 119961): the discovery of stars in the DB gap through the 
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Sloan Digital Sky Survey (SDSS; lEisenstein et al 
carbon-rich stars — the so-called Hot DQ stars ( iDufour et al.l 12007 



I), the discovery of a new class of 



of the DB stars in the ESO Supernova la Progenit or Survey (SPY 



far-ultraviolet observations of DB stars with FUSE ( Petitclerc et al.ll2005t iDesharnais et al. 



2008 



, the an alysis 



Voss et al. 2007), the 



20081 : iDufour et al.ll2010al ). the likelihood of a tmospheric contamin ation of the photospheres 
of DB stars through disk- fed accretion (e.g., IDufour et al.ll2010bl ). etc. Yet our knowledge 
of some of the funda r aental properties of DB stars, embodied in a set of questions posed 
by IBeauchamp et al.l (Il996l ). remains sketchy. These questions referred to the respective 
mass distributions for the DA and DB stars and to the possible existence of a bimodal 
mass distribution for DB degenerates, to the existence of statistical evidence for differences 
between the DB and DBA samples, as well as to the existence of atmospheric peculiarities in 
hot DB stars that might afford possible evidence for a recent episode of convective mixing. 



One particular problem with the analysis of DB white dwarfs is related to the effect of 
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invisible traces of hydrogen that can affect significantly the temperature s c ale of DB stars 
based on line profile fitting techniques. For instance, iBeauchamp et al.l (119991 see their 



Table 1) showed that small amounts of hydrogen, of the order of H/He ~ 10~^ — 10~^, can 
produce differences as large as 3000 K in Teg for the hottest DB stars. In the case of the 
Beauchamp et al. analysis, upper limits on the hydrogen abun dance were based only on 
the absence of H/3. To improve on this aspect of the analysis, IVoss et al.l ( 120071 ) secured 
high- resolution spectra in the region near Ha for a sample of ~ 60 DB stars in the SPY 



survey. A similar analysis had been carried out earlier by iHunter et al.l (120011 ) , but on a 
smaller sample of hot DB stars. With their high-resolution data, Voss et al. were able to 
secure more stringent limits on the hydrogen abundance. In their Table 1, they contrast the 
effective temperatures and surface gravities of DB stars without visible hydrogen obtained 
from two distinct grids of model atmospheres: one with a pure helium composition, and 
one with a trace of hydrogen, H/He = 10~^, the upper limit imposed by their data. As 
discussed in Voss et al., the average temper ature difference they me asure, around 230 K, is 
significantly smaller than that obtained by IBeauchamp et al.l ( 1l999l ) based on much higher 
upper limits on the hydrogen abundance. However, this evaluation is somewhat misleading 
since the upper limit imposed by the SPY data is strongly temperature dependent according 
to the results shown in Figure 9 of Voss et al. In fact, at Tcs ~ 24, 000 K, their detection 
limit of 300 mA for the Ha equivalent width indicates an upper limit of H/He ~ 10~^, or a 
factor of ten larger than the value used in their Table 1. Hence we feel that the influence 
of invisible amounts of hydrogen on the determination of the atmospheric parameters of DB 
stars still needs to be properly addressed. 

Another complication with the spectroscopic analysis of DB stars is related to the 
problei n of the convectiye effic iency in their atmospheres. While the spectroscopic anal- 
ysis of IBeauchamp et al.l (11999! ) relied on model a tmospheres calcul ated with the so-called 
ML2/a = 1.25 version of the mixing-length theory, IVoss et al.l (12007! ) adopted in their study 
a much less efficient a = 0. 6 version, i.e. the same convective efficiency as for DA white 



dwarfs (iBergeron et al.!ll995! ). This may represent an important source of discrepancy when 
comparing the atmospheric parameters from various analyses. While some authors argued 
that there is no reason to e xpect the mixing length parameterization to be different in DA 



and DB white dwarfs (e.g., ICastanheira et al.l 12006! ). we may as well argue that there is no 



reason to expect it to be identical. Hence we feel the problem of the convective efficiency in 
the atmospheres of DB stars needs to be reexamined as well. 

Of deeper astrophysical interest is the question of the possible evolutionary channels 
that could produce DB white dwarfs. The progenitors of DB white dwarfs are believed to 
be DA stars, transformed into helium-atmosphere white dwarfs as a result of the convective 
dilution of a thin radiative hydrogen layer (Mh ~ 10~^^ ^0) with the deeper convective 
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helium envelope (iFontaine fc Wesemaell 119871 ) . Evidence for this mechanism rests on the 
complete absence of helium-atmosphere white dwarfs in the Palomar-Green (PG) survey 
between T^s ~ 45, 000 K — where the coolest DO stars are found, and ~ 30, 000 K — 
where the hottest DB stars are located. This so-called D B ^ap has partially be en filled in 
by the discovery of extremely hot DB stars in the SDSS (lEisenstein et al.ll2006l ). although 
the number of DB stars in the gap remains a factor of 2.5 lower than what is expected based 
on the number of cooler DB stars in the SDSS. We are thus dealing with a DB deficiency, 
rather than with a real gap, in this particular range of effective temperature. How do these 
DB stars in the gap fit in our understanding of white dwarf evolution in general? 

All these questions remain current, and need to be addressed for a complete under- 
standing of the physical properties and evolutionary status of DB stars. With the hope of 
shedding some light on these issues, we present in this paper a complete and comprehensive 
spectroscopic analysis of over 100 DB stars. We first begin in Section 2 with a theoretical 
investigation of model atmospheres of DB stars to explore the observational and theoretical 
issues at stake. The observational data used in our analysis are then presented in Section 3. 
Our spectroscopic analysis follows in Section 4 and the results are interpreted in Section 5. 
We present our conclusions in Section 6. 



2. Theoretical Considerations 



2.1. Model Atmospheres 



Our model at mospheres and synthe t ic spe ctra are built from the model atmosphere 
code described in iTremblay fc Bergeron! (120091 ) and refere nces therein, in which we have 
incorporated the improved Stark profiles of neutral helium of lBeauchamp et al.l ( 119971 ). These 
detailed profiles for over twenty lines of neutral helium take into account the transition from 
quadratic to linear Stark broadening, the transition from impact to quasi-static regime for 
electrons, as well as forbidden c omponents. I n the cont ext of DB stars, these m odels are 
comparable to those described by iBeauchamp Jl995h and IBeauchamp et al.l (119961 ) , with the 
exception that at low temperatures (Tp ff < 10, 800 K ), we now use the free- free absorption 
coefficient of the negative helium ion of iJohnI (119941 ). For the hydrogen lines in DBA stars , 
we rely on the improved calculations for Stark broadening of iTremblay fc Bergeron! (120091 ). 
Our models are in LTE and include convective energy transport treated within the mixing- 
length theory. The theoretical spectra are calculated using the occupation formalism of 



Hummer fc Mihalad (119881 ) for both the hydrogen and helium populations and corresponding 
bound-bound, bound-free, and pseudo-continuum opacities. 
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An important issue related to the modeling of DB spectra is the inclusion of van der 
Waals broadening at low effective temperatures (Tp g < 15, OOP K), where Sta rk broadening is 
no longer the dominant source of line broadening. iBeauchamp et al.l ( 119961 . see their Figure 
6) show the effect of this broadening mechanism in a log g vs Tg^ diagram for DB stars with 
parameters determined from optical spectroscopy. The neglect of van der Waals broadening 
results in spuriously high log g values at low effective temperatures. In this series of synthetic 
spectrum calculations, we in c lude t he treatment of this broadening mechanism provided by 
Deridder fc Van RensbergenI (119761 ). 



Our model grid covers a range of effective temperature between Tcs = 11,000 K and 
40,000 K by steps of 1000 K, while the logg ranges from 7.0 to 9.0 by steps of 0.5 dex. 
In addition to pure helium models, we also calculated models with logH/He = —6.5 to 
—2.0 by steps of 0.5. In order to explore the effect of varying the convective efficiency on 
the predicted emergent fluxes, models have been calculated wi th the ML2 version of th e 
mixing- length theory, which corresponds to the prescription of iBohm &: Cassinellil (jl97ll ). 
with various values of a = i/H, i.e. the ratio of the mixing length to the pressure scale 
height, ranging from 0.75 to 1.75 by steps of 0.25. 



2.2. Exploration of the Parameter Space 



The determination of the atmospheric parameters, T^s, logg, and H/He, of DB stars 
represents an inherent difficulty in any analysis of their properties. While the atmospheric 
p arameters for DA sta rs derived from various studies agree generally well (see, e.g.. Figure 9 
of lLiebert et al.ll2005l ). a similar cor nparison fo r DB s tars is less than satisfactory, as shown 
for instance in Figures 3 and 4 of IVoss et al.l (120071 ). which comp are the values obtained 



from the SPY survey with those obtained from independent studies (IBeauchamp et al.lll999 
Friedrich et aPboOol : ICastanheira et al.l 120061 ). 



In addition to being rarer than their hydrogen-line DA counterparts, the hotter DB 
stars are characterized by an optical spectrum where the neutral helium transitions exhibit 
little sensitivity to effective temperature. As an illustration, we show in the left panel 
of Figure [1] the model spectra0 for various values of T^s (shown in blue), compared to a 
model template at T^s = 24,000 K (shown in red). There is a wide range of temperature, 
between Teg ~ 20, 000 K and 30,000 K, where the optical spectra are virtually identical. This 
underscores the need to secure well calibrated, high signal-to-noise spectroscopic observations 
to study these stars. Equivalently, we show in the left panel of Figure [2] the variation of the 



Unless specified, the models displayed here assume a value of a = 1.25. 
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equivalent width of He I A4471 as a function of effective temperature for various hydrogen 
abundances. In all cases, the variation of the equivalent width shows a wide plateau, about 
10,000 K wide, that inhibits the determination of accurate effective temperatures from the 
optical spectrum of DB stars in this particular range of temperature. Also shown in this 
figure are the measured equivalent widths in our sample of DB stars, arbitrarily located at 
Tes = 22, 000 K. Interestingly enough, the DB stars in our sample with the strongest lines 
have equivalent widths that are larger than predicted by any of our models shown in this 
panel. 

Luckily, we fare somewhat better with gravity indicators, as illustrated in the middle 
panel of Figure [TJ The large gravity- sensitivity of several segments containing neutral helium 
transitions (notably the 2^P — 6^D A3819 hne, the 4100-4200 A region that contains the 
2^P - 5^S A4121, 2^P - 6^ D A4144, and 2^P - 6 ^S A4169 lines, as well as the strong 
2^P — 5^D A4388 transition; IWickramasinghd Il979l ) can be used successfully t hrough the 



entire temperature range. We disagree in that respect with the statement of IVoss et al. 



( 120071 . see Section 3.1) that the dependence of line profiles on logg is smaller for the cooler 
DB stars. Interestingly, the variation of the He I A4471 equivalent width with Tefj is also 
strongly logg dependent, as can be observed in the middle panel of Figure |2J For low-gravity 
DB white dwarfs, the maximum is now strongly peaked near 22,000 K, while for high-gravity 
DB stars, the plateau becomes even wider, and the maximum equivalent width is reached 
at temperatures near 33,000 K. Again, variations of logg help very little in matching the 
maximum observed equivalent widths of the DB stars in our sample. 

Less well documented is the sensitivity of the optical spectrum of DB stars to the 
convective efficiency. Helium atmospheres of white dwarfs become convective below Tes ~ 
60, 000 K, but it is not until the temperature drops below ~ 30, 000 K that the helium 
convection zone becomes sufficiently irnporta nt, depending on the assumed convective ef- 
ficiency (see Figure 10 of iTassoul et al.lll990[ ). For completeness and future reference, we 
show in Figure |3] the extent of the helium convection zone in a 0.6 Mq DB white dwarf. 
Convective energy transport is generally treated in model atmosphere calculations within 
the mixing-length theory, and unless convection becomes adiabatic, the emergent model 
fluxes will depend on the assumed convective efficiency, parameterized in this case by the 
free paramet er a. This additio nal parameterization was first introduced in the context of 
DB st ars by iBeauchamd (19951 ) on the basis of work carried out in the hydrogen-line DA 
stars ( iBergeron et al.lll995l ). 



To illustrate this sensitivity, we compare in the right panel of Figure [T] model spectra 
calculated with a = 0.75 and 1.75. The effect is subtle, and only noticeable between Tes ~ 
18, 000 K and 28,000 K. At higher effective temperatures, the atmospheres of DB stars 
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become almost completely radiative (see Figure [3]), while at lower temperatures convection 
becomes adiabatic, and the convective efficiency no longer depends on the specific value of 
a. Despite the apparent lack of sensitivity to variations of a, the effects on the normalized 
line profiles are more pronounced. This is illustrated in the right panel of Figure [2], where 
the equivalent width of He I A4471 as a function of effective temperature is shown for model 
atmospheres calculated with various convective efficiencies. Once again, the line strength is 
sensitive to a for effective temperatures between 18,000 K and 28,000 K and, in contrast to 
the results shown in the previous panels, the largest predicted equivalent widths (achieved 
with a = 0.75) are larger than the largest measured values. This means that, were we to 
overestimate a in our models, all the stars would be shifted toward the maximum of the curve 
(i.e., near 30,000 K for a = 1.75) and form a clump. In contrast, were we to underestimate 
a, the stars we fit would be lumped on respective sides of the maximum of the curve (i.e., 
near 18,000 K and 30,000 K for a = 0.75), and gaps would develop in the distribution of 
stars with respect to effective temperature. This behavior can be used to constrain a, as we 
show below. 



We finally note that the results presented here are consistent with those of IVoss et al. 



(120071 ) who mention that the helium lines in their models, calculated with ML2/a = 0.6 (i.e. 
a value even smaller than the smallest value used in Figure [2]) , reach a maximum strength 
around T^s = 20, 000 K. 



3. Spectroscopic Observations 



Our sample of DB stars was selected from the online version of the Villanova White 
Dwarf CataloJ^. One of the first goals of this pr oject was to secure spectroscopic observations 
of all DB stars identified in the PG Survey (iGreen et al.l 119861) in order t o compare the 
luminosity function of DB stars with that obtained by iLiebert et al.l (120051 ) for DA stars. 
Otherwise, DB white dwarfs were routinely observed over the years and included in our 
sample. Due to observational constraints, we focused mainly on the brightest white dwarfs in 
the catalog, and as such, there is little overlap between our sample and the SDSS sample. The 
distribution of V mag nitudeflfor the 118 objects observed spectroscopically is presented in 
Figure m (we exclude here the 10 misclassified stars discussed further below). Several objects 
in this original sample are excluded from our analysis. Thre e are DA + DB double deg enerate 
systems and have been analyzed elsewhere: PG 1115+166 (IBergeron fc Liebertll2002l ). KUV 



^http://www. astronomy. villanova.edu/WDCatalog/index. html 

•^These can be photographic magnitudes in some cases, or even B magnitudes. 
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02196+2816 flLimoges et al.ll2009h . and KUV 14197+2514 f|Limoges fc Bergeronll2010h. Two 



are m agnetic DB stars: PG 0853+164 (jWesemael et al.ll200ll ) and Feige 7 (lAchilleos et al. 



19921 ). Finally, five white dwarfs do not show enough helium lines in their optical spectra to be 
analyzed quantitatively: LP 891-12 (0443-275.1), G102-39 (0551+123), CD 84 (0714+458), 
G227-5 (1727+560), and IDS 678A (1917-077). 

We are thus left with a total of 108 DB w hite dwarfs, including several DBA and DBZ 
stars (the analysis of iBeauchamp et al.l ( 119961 ) included close to 80 objects). The complete 
list of o bjects is provide d in Table 1. For comparison, a recent analysis based on the SPY 
sample (IVoss et al.ll2007l ). the largest completed up to now, includes 69 white dwarfs. Most of 
our objects were observed at high signal-to- noise ratio (S/N) in the blue region of the optical 
spectrum (3700-5200 A) with the Steward Observatory 2.3 m Bok Telescope equipped with 
the Boiler & Chivens spectrograph. The 4'.'5 slit together with the 600 line mm~^ grating 
blazed at 3568 A in first order provides a spectral coverage from about 3500 to 5250 A at a 
resolution of ~ 6 A FWHM. Additional blue spectra in the southern hemisphere were secured 
at Carnegie Observatories' 2.5 m Irenee du Pont Telescope at Las Campanas (Chile) with the 
Boiler & Chivens spectrograph. The 1'.'5 slit with the 600 line mm~^ grating blazed at 5000 
A provided a spectral coverage from about 3500 to 6600 A at a slightly better resolution of ~ 
3 A FWHM. Finally, the spectrum of GD 554 (2250+746) was obtained at the Mt. Hopkins 
6.5 m MMT telescope using the Blue Channel of the MMT Spectrograph. The 500 lines 
mm~^ grating and a I'.'O slit provided a spectral coverage from about 3400 to 6300 A with a 
resolution of ~ 4 A FWHM. 

Our blue optical spectra are displayed in Figure |5]in order of decreasing effective temper- 
ature. In addition to the numerous neutral helium lines that dominate the spectra of these 
stars, hydrogen lines (H/3) and the Ca ll H & K doublet are also visible in some objects. 
The spectrum of L7-44 (1708—871) suffers from residual flux calibration problem, while that 
of HE 0429—1651 (0429—168) is significantly contaminated by the pre sence of an M dwarf 
companion. Note that the spectra of the DB stars recently identified by iLepine et al.l (1201 if ) 
represent new observations at higher S/N. One of the most obvious qualities of our sample 
is its homogeneity, both in terms of wavelength coverage and signal-to-noise. While its size 
remains small compared to that of the SDSS sample, there is no comparison in terms of 
the quality of the spectra, as can be seen from Figure where we compare the distribution 
of S/N of our observations with that of the SDSS spectra from the Data Release 4. The 
majority of our spectra have S/N well above 50, while the SDSS spectra are strongly peaked 
between S/N ~ 5 and 20. Since the exposure time of a given SDSS spectrum is set by that 
of the entire plate, the corresponding S/N is necessarily a function of the magnitude of the 
star, and SDSS white dwarfs are much fainter than those studied here. 
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As part of our survey, we observed several DB white dwarf candidates that turned out 
to be lower surface gravity objects, mostly sdOB stars (with the glaring exception of one DA 
star). For completeness, we show these misclassified objects in Figure [71 Note that several 
of these have already been reclassified properly in the online version of the Villanova White 
Dwarf Catalog. 

Because coverage of the Ha line is essential for the determination of accurate atmo- 
spheric parameters for DB stars, we complemented our blue d ata with red spectra obtained 
from several distinct sources: the KPNO 4 m data used by iHunter et al.l (120011 ). spectra 
secured at the du Pont 2.5 r n telescope , spect ra extracted from the Data Release 4 of the 
SDSS, and SPY spectra from lVoss et al.l (120071 . only those with Ha detected), available from 
the ESO archive web site. We have also recently secured additional spectra at the Steward 
Observatory 2.3 m Bok Telescope and at the KPNO 4 m telescope. Our own spectroscopic 
observations (36 objects) as well as SDSS spectra (17 objects) are displayed in Figure [8] 
in order of right ascension; these have typically 3 to 5 A resolution (F WHM). SPY / UVE S 
spectra for 13 additional objects can be visualized as online material in IVoss et al.l ( 120071 ). 
To guide the eye, we also show in this figure the location of the Ha absorption feature. 
Hydrogen is detected in 50% of the obje cts shown here, a value in good agreement with the 
55% fraction of DBA stars estimated by lVoss et al.l (120071 ). This fraction drops to ~ 44% if 
we consider all stars in our sample, since 42 objects remain without Ha data. Note how the 
strength of Ha varies considerably from object to object, and how it is particularly strong 
in LP 497-114 (1311+129). 

Finally, in Secti on 4.4, we also rnake u se of /f/i? low-dispersion spectra i ssued from the re- 
ductio n procedure of lHolberg et al.l (120031 ) . which follows the prescription of lMassa fc Fitzpatrick 
(120001 ) for the correction of residual temporal and thermal effects and for the absolute flux 
calibration. A total of 34 DB stars in our sample have UV spectra suitable for analysis. 



4. Atmospheric Parameter Determination 
4.1. Fitting Technique 



Our fitting technique is similar to that outlined in iLiebert et al.l ( l2005l . and references 
therein) for the analysis of DA stars. The first step is to normalize the flux from individual 
predefined wavelength segments, in both observed and model spectra, to a continuum set 
to unity. The comparison with model spectra, which are convolved with the appropriate 
Gaussian instrumental profile, is then carried out in terms of this normalized spectrum only. 
The most sensitive aspect of this approach is to properly define the continuum of the observed 
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spectra. Here we rely on the procedure developed by lLiebert et al.l (120051 ) . where the fluxed 
spectrum is first fitted with model spectra, multiplied by a high order polynomial (up to A^), 
in order to account for any residual error in the flux calibration. T he nonlinear least-squares 
minimization method of Levenberg-Marquardt (jPress et al. 1986 ) is used throughout. Note 



that the values of Tes and logg derived in this first step are meaningless since too many 
fitting parameters are considered, and the model fit just serves as a smooth fitting function 
used to define the continuum of the observed spectrum. Normal points are then fixed at 
several wavelengths defined by this smooth function, and used to normalize the spectrum. 
This method turns out to be quite accurate when a glitch is present in the spectrum at the 
location where the continuum is set. It also provides a precise value of the line center, which 
can be corrected to the laboratory wavelength. An example of continuum fitting using this 
procedure is illustrated in the top panel of Figure M 

Once the spectrum is normalized to a continuum set to unity, as shown in the bottom 
panel of Figure [9], we use our grid of model spectra to determine Tes, ^ogg, and H/He in 
terms of the normalized spectrum only. Our three-dimensional minimization technique again 
relies on the nonlinear least-squares method of Levenberg-Marquardt, which is based on a 
steepest descent method. 



4.2. Determination of Tes, ^ogg, and H/He 

As shown in Figure [21 two solutions exist for a given DB spectrum, one on each side of 
the maximum strength of the neutral helium lines. In most cases, it is easy to distinguish 
these cool and hot solutions from an examination of our best fits, as displayed in the left 
panel of Figure [ID] for the cool DB star LP 475-242 (0437-M38). In addition, we also look 
at our spectroscopic solution in terms of absolute fiuxes (not shown here) by normalizing 
both the observed and predicted model spectra at 4250 A. In general, the slopes agree well, 
except in some cases where the object has been observed at high airmasses, in which case 
the slope of the observed spectrum might be distorted. 

However, when the star is close to the temperature where the helium lines reach their 
maximum strength — a temperature range that strongly depends on logg and on the as- 
sumed convective efficiency according to Figure [2] — the cool and hot solutions cannot be 
distinguished, as shown in the right panel of Figure UJJl for the hot DB star PG 1115+158. 
The comparison of the slopes of the energy distributions does not help either in this case. 
Fortunately, we discovered that this ambiguity could be resolved with the help of spectro- 
scopic observations at Ha, which add an additional constraint in our fitting procedure. Since 
our Ha data are independent of our blue observations, we use an iterative procedure where 
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we fit the blue spectrum with an assumed hydrogen abundance to obtain a first estimate 
of Tgff and log (7. The Ha spectra are then used to determine the hydrogen abundance, or 
upper limits on H/He, at these particular values of T^s and log (7. The procedure is then re- 
peated until the value of H/He has converged. An example of our solution for PG 1115+158 
is displayed in Figure [TTl The weak Ha absorption feature shown in the insert (from a 
SDSS spectrum, in this case) serves as an important constraint to determine the hydrogen 
abundance, which otherwise could only be inferred from the weak depression observed near 
H/3. Interestingly, our final solution for this object lies between the cool and hot solutions 
displayed in Figure [TOl 

The procedure described above works well when any hydrogen absorption feature (Ha 
or H/3) is visible in the optical spectrum. When no such feature is present, only upper limits 
on the hydrogen abundance can be determined. These limits depend on the signal-to-noise 
ratio of the observations, and on whether they are based on our blue spectra for H/3, or our 
red spectra for Ha. We show, in Figure [121 the limits on the hydrogen abundance imposed 
by the absence of Ha and H/3 features in our spectra (we adopt a detection limit of 200 
mA and 300 mA for the equivalent widths of Ha and H/ 3, respectiv e ly); th e results for Ha 



are qualitatively similar to those shown in Figure 9 of IVoss et al.l ( 120071 ). For high S/N 



spectra with no detectable hydrogen featur^, we set the hydrogen abundance in our fitting 
procedure at the value given by these upper limits at the appropriate temperature. In cases 
where the spectrum is noisier, our fitting procedure may find an upper limit on the hydrogen 
abundance that is larger than those of Figure [T2i 



4.3. DBZ White Dwarfs 

Four DBZ stars in our sample have sufficiently strong calcium lines that the inclusion of 
the Ca II H & K doublet is required to fit them properly. These are GD 40 (0300-013), CBS 
78 (0838+375), KUV 15519+1730 (1551+175), and G241-6 (2222+683). We take a simple 
approach here of including only calcium in our equation of state and opacity calculations. 
A more detailed analysis of these DBZ stars, including additional heavy elements, will be 
presented elsewhere. A smaller grid of models between T^s = 12, 000 K and 17,000 K has 
been calculated for this purpose, with calcium abundances of log Ca/He = —7.5 to —6.0 
by steps of 0.5, and the same logg and H/He values as above. For each star we adopt the 



^W e use in particular the fact that no detectable Ha absorption feature has been reported bv lVoss et al. 
(l2007i see their Table 2) in 0308-565, 0429-168, 0845-188, 0900+142, 1046-017, 1144-084, 1336+123, 
2129+000, and 2144-079. 
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calcium abundance that best reproduces the calcium doublet. Our best fits for these 4 DBZ 
stars are presented in Figure [131 



4.4. Convective Efficiency in DB White Dwarfs 



T he problem of the cony ective efficiency in the atmosphere of DA stars has been tack- 
led by iBergeron et al.l ( 119951 ) who used optical spectroscopic observations combined with 
UV energy distributions to show that the so-called ML2/a = 0.6 parametrization of the 
mixing-length theory provides the best internal consistency between optical and UV effective 
temp eratures, trigonorn e tric p arallaxes, V magnitudes, and gravitational redshift^. Simi- 
larly, iBeauchamp et al.l ( 1l996l . see their Figure 5) showed a comparison of optical and UV 
temperatures of DB white dwarfs using model atmospheres calculated with the MLl, ML2, 
and MLS versions of the mixing-length theory; the authors only concluded then that ML2 
and ML3 provided a mo re satisfactory correlation than MLl. The spectroscopic analysis of 
Beauchamp et al.l ( 1l999l ). however, relied on model atmospheres calculated with a convective 
efficiency of ML2/a = 1.25, although details of this particular choice were not provided in 



that s tudy. The spectroscopic analysis of the DB stars in the SPY survey by IVoss et al. 
( l2007l ). on the other hand, relied on models calculated with ML2/a = 0.6, i.e. the same con 



vective efficiency as for DA models. As shown in Figure [H the strength of the helium lines 
in DB stars are significantly affected by the assumed parameterization of the mixing-length 
theory used in the model atmosphere calculations, in particular between Teg ~ 18, 000 K 
and 30,000 K. It is thus of utmost importance to calibrate the convective efficiency in DB 
stars as accurately as possible. 



The log g versus T^g distribution of all DB stars in our sample is displayed in Figure 
for three versions of the mixing-length theory, namely ML2/q; = 0.75, 1.25, and 1.75. As 
expected, only objects with 18,000 K < T^s < 30, 000 K are affected by this parameterization. 
As discussed at the end of Section 2.2, if the convective efficiency is underestimated, the 
helium lines are predicted too strong near the maximum equivalent widths, and DB stars 
will have a tendency to be lumped on either side of this maximum. This is precisely what is 
observed in Figure [H] for a = 0.75 where the max imum occurs near 23,000 K according to 
Figure [2] (see also Figure 3 of iBergeron et al.lll995l for DA stars near the ZZ Ceti instability 
strip). On the other hand, if the convective efficiency is overestimated, helium lines are 
predicted too weak near the maximum, and DB stars in this case will tend to accumulate at 



'"'The more recent analysis by Tremblav fc Bergeron ( 20091 ) . based on improved Stark profiles, suggests . 
more efficient value of a = 0.8. 
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this location. Again, this is what is observed in Figure [TUfor a = 1.75 where DB stars form a 
clump at Tcfj ~ 25, 000 K, even leaving a void near 20,000 K (see a similar result for DA stars 
in Bergeron et al. 1995). The results with a = 1.25, however, provide a smoother distribution 
and a uniform increase of the number of DB stars as a function of effective temperature, 
which is exactly what is expected in terms of the white dwarf luminosity function. Based 
on these results alone, a value around a = 1.25 seems entirely reasonable, even though we 
cannot easily rule out values between a ~ 1.0 and 1.5 at this stage. 



Bergeron et al.l (119951 ) went a step further to calibrate the convective efficiency in DA 
stars, by comparing optical temperatures determined from spectroscopy with those obtained 
from UV energy distributions. Such a comparison is displayed in Figure [15] for two versions 
of the convective efficiency, namely ML2/a = 1.25 and 1.75. Also shown in the bottom panel 
is a comparison of our UV tem peratures, derived f r om fit s to spectra with ML2/a = 0.6 



models, with those obtained by lCastanheira et al.l (120061 ) based on similar assumptions; note 



the almost perfect agreement here. In all UV fits, we simply assume log g = 8 and pure helium 
compositions for all objects. We first begin by examining the results for a = 1.25. Although 
there is a good overall agreement between optical and UV temperatures, there are also several 
discrepant results at the hot end of the sequence. The two obj ects with Tppt ^ 30, 00 K 



correspond to PG 0112+104, recently analyzed in great detail by iDufour et al.l (l2010al ). and 
PG 1654+160. Although the /f/i? spectrum of the first object is of good quality, the spectrum 
of the second object suffers from a significant discontinuity between the SWP and LWP 
images. It is also possible that the UV energy distributions may become increasingly less 
sensitive to temperature at the hot of the DB sequence, while interstellar extinction may have 
to be taken int o account for t h e intrin sically brighter, and more distant, DB stars. Indeed, for 



PG 0112+104. iDufour et al.l (l2010al ) found the energy distribution extending into the FUV 
region (~ 925A) globally consistent with the ~ 31, 000 K effective temperature derived from 
the optical spectroscopy provided that a small amount of reddening were included in the fits 
{E{B -V) = 0.015). The most important outlier in Figure US] is LP 497-114 (1311+129) 
with a ~ 8000 K temperature difference; this peculiar white dwarf is further discussed in 
Section 5.4. 

More worrisome in the top panel of Figure [TS] is the systematic offset between both 
temperature scales (with Tuy > Topt) for 15,000 K < Tyy < 22, 000 K. Note that a more 
efficient prescription of the mixing-length theory, shown in the middle panel of Figure [T5| 
improves the agreement between both temperature scales near Tyy ~ 22, 000 K, but the 
agreement gets worse for hot DB stars, and the systematic shift still remains for 15,000 K 
< Tuy < 20, 000 K. Hence, it is not clear that increasing the convective efficiency would 
improve the overall agreement. 
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Since UV temperatures may depend on surface gravity, hydrogen abundance, and con- 
vective efficiency, we explore in Figure [16] the effects of variations of these parameters on our 
temperature determinations. All in all, we see that UV temperatures are fairly independent 
of our initial assumption of log g = 8 and pure helium compositions in our fit^. The same 
conclusion applies to our particular choice of a. We thus conclude that the systematic shifts 
observed in Figure [15] with our ML2/a = 1.25 models cannot be explained in terms of inad- 
equate assumptions in our fits to UV data. What the overall results suggest, instead, is that 
theoretical improvements probably need to be made at the level of the model atmospheres, 
perhaps even with the treatment of convective energy transport (see below). 



4.5. Error Estimation 



The internal uncertainties of the atmospheric parameters — T^s, log (7, and H/He — 
can be obtained directly from the covariance matrix of our ffiting method. These depend 
mostly on the signal-to-noise ratio of the observations and on the sensitivity of the models 
to each parameter. For DA white dwarfs, these internal errors can become negligibly small 



with sufficiently high sig nal-to-noise spectroscopic observations (S/N > 50; [Bergeron et al. 



1992 . [Liebert et al.ll2005l ). Since our DB spectra are rather homogeneous in terms of S/N (see 
Figures [3 and [HD and of sufficiently high quality (S/N > 50), the internal errors in our analysis 
will be mostly dominated by the sensitivity of the models to each fitted parameter. The 
theoretical spectra displayed in Figure [1] indicate that the internal errors will be relatively 
small for log (7 (similar conclusions apply to H/He), but can be relatively large for T^g, 
particularly in the temperature range where the strength of the helium lines reaches its 
maximum. 

The true error budget, however, must also take into account external uncertainties 
originating mostly from flux cali bration, which can b e estimated from multiple observations 



of the same star, as performed by [Liebert et al.[ (|2005l . Figure 8) for DA stars, or by [Voss et al. 



(|2007[ . Figures 1 and 2) for DB stars. Even though we have not reobserved DB stars for that 
specific purpose (except in a few cases), we still managed to secure repeated observations of 
28 DB stars in our sample, either because of flux calibration problems (e.g., objects observed 
at high airmass), or insufficient S/N. Hence the external errors determined below will tend 
to be overestimated, if anything. 



^In light of these results, we are surprised that Castanheira et al. ( 20061 ) succeeded in constraining logg by 
using only lUE data (se e their Tabl e 1); n ote that their log g determinations were omitted in the comparison 
displayed in Figure 4 of [Voss et al.[ ([20071 ). 
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We compare in Figure [T7I the T^s and log (7 measurements for the 28 DB stars in our 
sample with mutliple spectra. As expected, the external error on Teg are particularly large 
between ~ 20, 000 K and 25,000 K, where the helium lines reach their maximum strength. 
Otherwise, outside this temperature range, the temperature estimates are in excellent agree- 
ment. The spread in log (7 measurements is more significant, however, but some of the most 
extreme outliers are white dwarfs with a first spe c trosco pic observation at very low S/N; these 
include the DB stars reported by iLepine et al.l (120111 ) for instance. If we take the average 
uncertainties of both parameters, we obtain {ATes/T^s) = 4.6% and (Alogf?) = 0.069, 
but if we remove the obvious outliers, these numbers drop to {ATcs/Tcs) = 2.3% and 
(A log (7) = 0.052. Given that these are probably overestimated, as discussed in the previous 
paragraph, we adopt these values as conservative estimates of the external uncertainties of 
our atmospheric parameters. Note that these results are comparable to those obtained by 
Voss et al.l J2007h . (AT^s/T^s) = 2.03% and (Alog^) = 0.058, based on multiple observa- 



tions of 24 objects. 



4.6. Adopted Atmospheric Parameters 

The atmospheric parameters for the 108 DB and DBA stars in our sample are reported 
in Table 1; the values in parentheses for T^g and logg represent the combined (in quadra- 
ture) internal and external errors of the fitting technique, while only the internal error is 
available for H/He. For DB stars without detectable hydrogen, upper limits on the hydrogen 
abundance are given in Table 1. The stellar mass (M) and white dwarf coolin g age (log r) of 



each s tar are obtained from evolutionary models similar to those described in lFontaine et al. 



fl200lh but with C/0 cores, g(He) = log Mue/M^ = 10"^ and g(H) = lO^^", which are 
representative of helium-atmosphere white dwarfs. The absolute vis ual magnitude (My ^ 



and lu minosity (L) are determined with the improved calibration from iHolberg fc Bergeron 



(120061) . defined with the Hubble Space Telescope absolute fiux scale of Vega. The V magni- 
tudeCl are taken from the Villanova White Dwarf Catalog, which combined with My yield 
the distance D. Finally, for all PG stars in the complete sample, we also provide the 1/Knax 
weighting (pc~^) used in the calculation of the luminosity function, where Vmax represents 
the volume defined by the maximum distance at which a gi ven object would st ill appear in 



the sample given the magnitude limit of the PG survey (see lLiebert et al.ll2005l for details). 



Sample fits for both DB and DBA stars covering the full temperature range of our sample 



^Note that these sometimes represent photographic magnitudes, Stromgren y magnitudes, or even B 
magnitudes. 
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are displayed in Figure [181 The left panels show the blue portion of our spectroscopic fits, 
while the right panels show the corresponding region near Ha. The importance of the red 
coverage is clearly illustrated, with Ha being barely detected in several DBA stars, while 
H/3 remains spectroscopically invisible. In other objects, only an upper limit on H/He could 
be set, based on the absence of Ha. We also show in this figure the DBA star with the 
largest hydrogen abundance measured in our sample (LP 497-114; 1311+129). A trend for 
the coolest objects to have larger than average surface gravities is already apparent. 



5. Selected Results 
5.1. Comparison with Results from SPY 

As mentioned in the Introduction, th e spectroscopic analysis of ~ 70 DB and DBA 



stars in the SPY survey (jVoss et al.l 120071 ) represents the largest set of data against which 
our atmospheric parameter determinations can be compared. We have 44 white dwarfs in 
common with the SPY sample, 22 of which are DBA stars. The comparison of effective 
temperatures and surface gravities is displayed in Figure 



Voss et al.l (120071 ) also included van der Waals broadening in their models, in a simplified 
form. Even though they found a much better agreement with the observed line profiles and 
strengths, their fits were not completely satisfactory according to the authors; they also men- 
tion that a fit in which log g is actually allowed to vary converges to very high values. Hence 
a value of log g = 8 was simply assumed for the coolest DB stars in their sample (shown as 
open circles in Figure [T9l) . The comparison of logg values is thus meaningless for these stars, 
although we note that the corresponding temperatures are in excellent agreement, except 
for the two coolest white dwarfs for which our values of log ~ 9 are l ikely overestimated . 



Worth mentioning in this context are the results shown in Figure 5 of iKepler et al.l (120071 ) 
for the DB stars identified in the Data Release 4 of the SDSS, where the masses gradually 
increase to very large values (M > 1.0 Mq) at low effective temperatures. Whether this 
increase is due to an inappropriate treatment of van der Waals broadening, or to the neglect 
of this broadening mechanism altogether, is not discussed in their analysis. 

For Tes < 19, 000 K, the T^s determinations displayed in Figure [19] are in good agree- 
ment, although our logg values in this temperature range appear systematically larger than 
those of Voss et al. by 0.15 dex, on average. The logg values agree much better at higher 
temperatures, however, with the exception of PG 1115+158 (labeled in the figure) for which 
we get a surface gravity 0.4 dex higher than the value derived by Voss et al. (logg = 7.52). 
This white dwarf is of particular interest because of its low inferred mass of only 0.385 Mq, 
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and also because it is a member of the pulsating V777 Her class. iBeauchamp et al.l (Il996[ ) 
found no low-mass DB stars in their sample and argued that evolutionary scenarios that 
could produce such low-mass degenerates could simply not form DB white dwarfs. Both 
spectroscopic solutions for PG 1115-1-158 are compared in Figures [TT] and |2UJ While the 
effective temperatures agree well within the uncertainties, the low surface gravity obtained 
by Voss et al. predicts helium lines that are much shallower than observed. Hence we believe 
that our solution at log (7 = 7.91 (or M = 0.56 Mq) is more appropriate for this star. 

The differences in temperature for T^g > 19, 000 K are significantly more important, 
partly due to the use of different versions of the mixing-length theory (ML2/a = 0.6 versus 
a = 1.25), but also because of the intrinsic difficulty of measuring Tcs precisely in this par- 
ticular temperature range. The most extreme case here is for KUV 03493-f 0131 (0349-1-015; 
labeled in Figure [TI?]) for which we obtain an effective temperature ~ 6000 K higher than 
Voss et al. Both spectroscopic solutions for this object are displayed in Figures [18] and [201 
Note that if we rely only on the blue spectrum for this object, we find a cooler solution 
at Tcs = 22, 760 K, still significantly hotter than Voss et al. We find that our solution at 
Tgff = 24, 860 K reproduces the widths and overall strengths of the neutral helium lines, as 
well of the observed slope of the spectral energy distribution (not shown here), much better 
than a solution at 18,740 K. 



5.2. Mass Distribution 



Beauchamp et al.l (119961 ) showed that the mass distribution of DB stars was relatively 
narrow, with the mass of 73% of their 41 DB stars above Teg = 15, 000 K (excluding DBA 
stars) lying between 0.5 and 0.6 Mq, for an average mass of (Mdb) = 0.567 Mq. The 13 
DBA stars were found at slightly larger masses, with an average mass of (Mdba) = 0.642 
Mq, suggesting that more massive white dwarfs may have a tendency to show hydrogen, 
perhaps through an increased hydrogen accretion from the interstellar medium, or through 
a decreased dilution of t he accreted hydrogen throughout a thinner helium convection zone 
( IBeauchamp et al.lll996l ). One significant distinction between the mass distributions of DB 
and DA stars, noted by Beauchamp et al., was the almost complete absence of the low- 
and high-mass tails in the mass distribution of DB stars, suggesting that the formation 
of double degenerates and mergers, which are often invoked to explain respectively these 
low- and h i gh-ma ss tails for DA stars, are not producing DB stars. The picture drawn by 
Voss et al.l (120071 ) using ~ 50 DB stars fror a the SPY survey differs significantly from the 
one described above. While the results of [Beauchamp et al.l (Il996[ ) indicate that 25% of 
all DB white dwarfs are DBA stars, the high-resolution spectroscopic observations at Ha 



- 18 - 



of the SPY survey revealed a much larger fraction of 55% of DBA stars in their sample, 
as discussed above. The resulting mean mass of each subsample, (Mdb) = 0.584 Mq and 
(Mdba) = 0.607 Mq, are now in much closer agreement, suggesting that both populations 
have a common origin. Furthermore, Voss et al. found a significantly larger mass dispersion, 
with several low-mass (M < 0.5 M©) and high-mass white dwarfs in their sample, in line 
with the results for DA stars. 

The distribution of mass as a function of effective temperature for all 108 DB and 
DBA stars in our sample is displayed in Figure |2T] Our final sample comprises 47 DBA 
stars, or 44% of the DB white dwarf population. Based on our discussion above, it is 
clear that this fraction represents only a lower limit since we are lacking the red spectral 
coverage for many of the DB stars in our sample, the cool o nes in particu l ar. In contrast 



to the mass distribution of DB stars in the SDSS obtained by iKepler et al.l ( 120071 . see their 
Figure 5), our mass distribution does not show this overwhelming increase in mass at low 
effective temperature, even though we do find several stars with masses in excess of 1 Mq. 
These are (below Teg = 15,000 K) 0249+346, 1419+351, 1542-275, 2058+342, 2147+280, 
and 2316—273, all of which are shown in the last panel of Figure Note that there are 
additional white dwarfs in the same range of temperature that have more normal masses 
(M < 0.8 M0), for instance 0517+771, 1056+345, 1129+373, 1542+244, also shown in this 
last panel, with the exception that these have well defined helium absorption features — the 
log 5f-sensitive He I A3819 line in particular — as opposed to the very massive ones. Instead 
of invoking a problem with the model spectra, our results suggest that we have probably 
reached the limit of the spectroscopic technique for these objects. They are probably cooler, 
more normal mass DB stars, with extremely weak absorption features (see further evidence 
below). 

Even if we ignore the most massive objects in Figure |2T1 we still note a trend for DB 
stars in the range 13, 000 K < T^q < 18, 000 K to have larger masses, up to ~ 0.9 Mq 
i n some cases, a problem usually attributed to the treatment of van der Waals broadening 



( iBeauchamp et al.lll996l ). What is more intriguing, however, is that white dwarfs with normal 
masses, M ~ 0.6 Mq, are also found in the same range of temperature as these massive stars, 
which implies that the large dispersion in mass that occurs at these temp e rature s might be 



real after all. This point has already been made by iLimoges &: Bergeron! (|2010[ ) who show 
in their Figure 7 the spectroscopic fits for two DB stars with nearly identical temperatures 
(Teff ~ 15, 000 K), but with log (7 values that differ by 0.33 dex (logg = 8.03 and 8.36). The 
comparison reveals that the He I A3819 and A4388 lines, which are the most gravity sensitive 
in this range of temperature, differ markedly in both stars, indicating that the high logg 
value measured in one of these stars is probably real, and not a simple artifact produced by 
the models. 
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The problem thus rests on our abihty to confirm independently the atmospheric param- 
eters obtained here within the current theoretical framework. To do so we rely on indepen- 
dent distance estimates obtaine d from trigonometric p arallax measurements taken mostly 



from the Yale parallax catalog (Ivan Altena et al. 



I994J, he reafter YPC) and from the new 



Hipparcos-h&sed parallaxes from Gould fc Chanam^ ( 2004 ). We found good trigonometric 



parallaxes for 11 DB white dwarfs in our sample. The parallax and corresponding distance 
for each object are reported in Table 2 together with our atmospheric parameter solution, 
including the spectroscopic distance. The comparison of spectroscopic distances and those 
obtained from trigonometric parallaxes is displayed in Figure [221 The agreement is surpris- 
ingly good, especially given the fact that white dwarfs with the closest match are found in 
the temperature range where the mass dispersion is large in Figure [21] (Teg ~ 15, 000 K; see 
Table 2). We are thus fairly confident that our surface gravity and mass values measured 
spectroscopically are fairly accurate. 

We also observe some discrepant distance estimates, labeled in Figure [221 For Feige 4 
(0017+136; label 1), we obtain a spectroscopic distance twice as large as that inferred from 
the trigonometric parallax. Our fit to this DB star is excellent, and we find no easy way 
to reconcile the two distance estimates. GD 358 (1645+325; label 2) is a pulsating white 
dwarf, the only hot DB star in Table 2. Again we do not have an easy explanation for the 
discrepancy. G188-27 (2147+280; label 3) has a spectroscopic distance significantly smaller 
than that inferred from the parallax. In this case, however, the spectroscopic log g value of 
8.85 (or M = 1.12 Mq) is certainly overestimated, and this object corresponds to one of 
the (almost) featureless DB stars discussed earlier. A value of log (7 = 8.2 would actually 
reconcile both distance estimates perfectly. This corresponds to a mass of ~ 0.7 Mq, i.e. the 
average mass of the bulk of DB white dwarfs near 15,000 K. 

The mass distribution of all DB and DBA stars in our sample, regardless of their tem- 
perature, is displayed in Figure [231 If we exclude from this distribution the most massive 
objects near ~ 1.2 Mq, which all correspond to the almost featureless cool DB stars dis- 
cussed above, we find a mean mass for our sample of (M) = 0.671 Mq with a standard 
deviation of only ctm = 0.085 Mq. In contrast to our preliminary conclusion presented in 



Beauchamp et al.[ (Il996l ) , we now find that there is no significant mass difference between the 



DB and the DBA samp les — (Mnp ) = 0. 657 Mq and (Mdba) = 0.688 Mq — in agreement 



with the conclusions of [Voss et al.l (120071 ). This result was already apparent from a quick 
examination of Figure [211 Note that the hottest white dwarfs in this figure are all of the DB 
type, and that their masses are slightly below 0.6 Mq, thus contributing to an average mass 
for the DB white dwarfs that is only slightly lower than the mean for DBA stars. 



We also compare, in Figure [231 the mass distribution of DB stars with that of DA stars 
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from the P G survey, reanalyzed with ou r improved models for hydrogen-atmosphere white 
dwarfs 
(Mda) 



,see 



Tremblay fc Bergeron! l201ll and references therein), and for which we obtain 
0.631 Mq and cta/ = 0.133 Mq. While the peak of the mass distribution of DB 
stars agrees with that of DA stars, the former distribution is shifted towards higher masses, 
resulting in an average mass ~ 0.03 Mq higher. The peak value for the DB stars , between 
0.60 and 0.65 Mq, also agrees with the peak value determined by lVoss et al.l (120071 see their 
Figure 8), although their mean mass of (M) = 0.596 Mq is much smaller than ours. This 
result is consistent with our logg determinations being larger than their values by about 
0.15 dex in the temperature range where the bulk of DB white dwarfs is found (see Figure 

m. 



Also of interest is the absence of low-mass DB stars in our sample, in contrast with 
low-mass DA stars, which are present in la rge numbers in all surveys, including the PG 



survey shown in Figure [231 Voss et al. ( 2007 , see their Figure 8) identified three object^ in 



the SPY sample with M < 0.5 Mq, the lowest mass DB white dwarf being PG 1115+158, 
with a mass of only 0.385 Mq (or logg = 7.52). This object has already been discussed 
in Section 5.1, and we find instead a significantly larger value of M = 0.56 Mq for the 
same star. Similarly, Voss et al. report a mass of M = 0.481 Mq (Teg = 16, 904 K) for 



H E 0423-1434, whi 



of iBeauchamp et al 



e we fi nd 0.64 Mq [Tes = 16,900 K). Hence we reaffirm the conclusion 
(I19961) that low-mass D B stars are rare, or absent, a conclusion also 



reached by [Bergeron. Leggett. &: RuizI (120011 ). These authors found, in their photometric 
analysis of 152 white dwarfs with trigonometric parallax measurements, that all low-mass 
degenerates probably possess hydrogen-rich atmospheres. As discussed in that study, since 
common envelope evolution is required to produce white dwarfs with M < 0.5 Mq — the 
Galaxy being too young to have produced them from single star evolution — we must 
conclude that this particular evolutionary channel does not produce helium-rich atmosphere 
white dwarfs. This could be the case because the objects which go through this close-binary 
phase end up with hydrogen layers too massive to allow the DA to DB conversion near 
Teff ~ 30, 000 K, or below. 



5.3. Luminosity Function 



Liebert et al.l (120051 ) presented spectrophotometric observations of a complete sample of 



348 DA white dwarfs identified in the Palomar Green Survey and obtained robust values of 
Teff, logs', masses, radn, and cooling ages for all stars in their sample using the spectroscopic 



Note that we could only find two white dwarfs with M < 0.5 Mq in their Tables 1 and 2. 
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technique. They also calculated the luminosity function of the sample, weighted by l/\4iax, 
where Knax represents the volume defined by the maximum distance at which a given object 
would still appear in the sample given the magnitude limit of the PG survey. An overall 
formation rate of white dwarfs in the local Galactic disk of 1 it 0.25 x 10~^^ pc^'^ yr^^ was 
also reported. 

Since we observed all DB stars in the PG survey, we can calculate their luminosity 
function in the same fashion as for DA stars using the 1/Knax values given in Table l|j. 
Our results are reported in Table 3 and shown in Figure [2H Note that we rely here on 



bolometric magnitudes rather than absolute visual magnitudes — as iLiebert et al.l (120051 ) 
did (see their Figure 10) — since DA and DB stars have different My values at a given 
effective temperature and surface gravity. Also shown in Figure [21] is the contribution from 
the DA stars in the PG survey, again reevaluated using our updated grid of DA model 
spectra. In the bolometric magnitude range over which DB stars are detected (roughly 
A^boi = 7— 11), the DB/(DA+DB) ratio is of the order of 0.2. We note, however, a large and 
significant increase in the value of log in the range from Mboi = 9 to 10, which corresponds 
to an effective temperature near 20,000 K. Thus, the DB/DA ratio is lower than the nominal 
value of 1 out of 4 white dwarfs at high effective temperatures, but increases sharply for stars 
below 20,000 K. This, it turns out, is also the temperature at which the bottom of the helium 
convection zone sinks sharply into the envelope along the cooling sequence (see Figure [3]). It 
is also the temperature below which DBA stars start to appear in large numbers in Figure 
[m and also where the average mass of all DB stars increases by ~ 0.1 Mq. 

The total space density of DB stars in the PG survey can be obtained by summing the 
values in Table 3 over all magnitudes bins. We obtain a space density of 5.15 x 10~^ DB 
star per pc^. Hence we expect in a volume defined within 20 pc from the Sun around tw o 



genuine DB stars, yet none have been identified in our sample (see also lHolberg et al.ll2008l ). 



5.4. Hydrogen Abundance and Mass Fraction in DB White Dwarfs 

About 44% of all DB stars in our sample show traces of hydrogen, or 50% if we consider 



only t he objects for which we have spectroscopic data covering the Ha region. IVoss et al. 



( 120071 ) obtained an even higher fraction of 55% based on the SPY data. Hence the DBA 
phenomenon is quite common among DB stars. The hydrogen abundances as a function of 
effective temperature for all DBA stars in our sample are displayed in Figure |25l Also shown 



^We ignore here the magnetic DB star PG 0853+163 as well as the DB component of the double degenerate 
binary PG 1115+166. 
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are the upper limits on the hydrogen abundance for DB stars, as determined from the absence 
of Ha or H/3. In general, these DB stars are ahgned on the observational limits reproduced 
here from Figure [T21 but some objects have noisier data and these limits are simply not 
reached. The first striking result is that almost all white dwarfs above T^s ~ 20, 000 K 
are DB stars, with the exception of 3 objects (KUV 05134+2605, PG 1115+158, and PG 
1456+103). 

There are two possible origins for hydrogen in these helium-dominated atmospheres. 
First, hydrogen could be residual, originating from the thin hydrogen atmosphere of the DA 
progenitor, convectively diluted into the more massive heliur n envelope during the DA to 



DB transition that occurred at higher effective temperatures. iMacDonald &: Venned (Il99ll ) 
showed that the exact temperature at which this process occurs depends on the thickness 
of the hydrogen layer (see their Table 1). For their SI models, which correspond to the 
Schwarzschild criterion with a value of a = 1, a DA star with hydrogen layer masses of 
Mh = 10-^^ 10-^^ and lO'^^ would respectively mix at temperatures of Teg = 27, 400 K, 
17,900 K, and 11,700 K. The main effect of adding more hydrogen on top of the helium con- 
vection zone is to delay the time it takes for helium to develop a sufficiently deep convection 
zone, thus reducing the temperature at which mixing occurs. Upon mixing, it is then as- 
sumed that the thin hydrogen atmosphere will be thoroughly mixed with the more massive 
helium convection zone, and that the resulting photospheric hydrogen abundance will sim- 
ply correspond to the ratio of the total hydrogen mass Mh to that of the helium convection 
zone, Mne-conv One problem immediately arises when trying to account for the hot DB 
stars above T^s ~ 20, 000 K in Figure [23 Indeed, for convective mixing to occur at these 
temperatures, the required hydrogen layer mass has to be of the order of 10~^^ Mq according 
to MacDonald & Vennes. This mass, it turns out, is comparable to the mass of the helium 
convection zone in this temperature range (see Figure [3]). The complete dilution of hydrogen 
should thus lead to a H/He ratio ~ 1 instead of the observed limits of H/He ~ 10"'^ — 10~^. 

We can explore this problem more quantitatively by converting the hydrogen abundances 
from Figure [25] into total hydrogen masses by assumin g that hydrogen is i ndee d homoge- 
neou sly mixed in the helium convection zone (see also 



Dufour et al 



2007a and Voss et al. 



20071 ). For simplicity, we assume evolutionary models at 0.6 Mq (as in Figure [3]) and we also 
ignore the feedback effect of the presence of hydrogen on the depth of the mixed H/He con- 
vection zon The total hydrogen mass as a function of effective temperature for all DBA 
stars in our sample is displayed in Figure [261 Again, the values for DB stars represent only 



^"The values of A/h inferred here differ somewhat from the values obtained bv iDufour et al.l (j2007al see 
Figure 13) since in their calculations, the depth of the helium convection zone w as taken from Table 1 of 



Dupuis et al. ( IQQSh . while we rely here on the more recent evolutionary models of Fontaine et al. ( 2001 ). 
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upper limits. Above Tcs = 22, 000 K, these upper limits imply Mh < 10"^^ Mq. Such small 
hydrogen layer masses are simply unable to maintain a hydrogen-rich atmosphere at higher 
temperatures. The only way to account for the existence of these hot DB stars in our sample 
is thus to conclude that they must have maintained a helium- dominated atmosphere through 
their entire life history. This interpretation is consi stent with the existenc e of hot DB stars 
in the DB gap, i.e. 30,000 K < Teg < 45,000 K jEisenstein et al.lbooeh . Even the three 
DBA stars near 24,000 K have hydrogen layer masses (Mh ~ 10~^^ Mcp) too small to have 
been D A stars in the past according to the results shown Figure 1 of iMacDonald fc Vennes 
( ll99ll ). However, our estimates are sufficiently approximate that larger masses of the order 
of Mh ~ 10^^^ Mq can probably be accommodated in these objects, in which case they 
would have undergone the D A to DB transition near ~ 27, 000 K according to Table 1 of 
MacDonald fc VennesI Jl99lh . 



We note that a similar evolutionary channel is required to explain the existence of the 
Hot DQ stars, w hose atmospheres a re dominated by carbon, with only traces of helium 
and no hydrogen (IDufour et al.ll2008l ). Since the hottest DQ stars known have temperatures 
around 24,000 K, their immediate progenitors must necessarily be DB stars with no hydrogen. 
Hence it is tantalizing to suggest that some of the hot DB stars observed here, as well as the 
hot DB stars in the gap identified in the SDSS, are the immediate progenitors of the Hot DQ 
stars. If this interpretation is correct, what fraction of DB stars between Tgff ~ 30, 000 K 
and 24,000 K are progenitors of Hot DQ stars remains to be determined. 

The situation is even more complicated below T^g ~ 20, 000 K where the inferred hy- 
dr ogen layer masses in th e DBA stars are in the range 10~^^ < Mh/M© < 10~^°. According 
to iMacDonald fc VennesI (119911 ). DA stars with such large amounts of hydrogen would only 



mix at temperatures below 12,000 K. In other words, we observe much more hydrogen in 
those stars than would be expected on the basis of the complete mixing of the hydrogen 
layer in the underlying helium convection zone. Hence, the hydrogen abundances in these 
DBA stars are too high to have a residual origin, and external sources of hydrogen must 
be invoked, either from the interstellar medium or from other bodies such as comets, dis- 



rupted asteroids, small plan ets, etc., a conclusion also reached by lVoss et al.l (120071 . see also 
MacDonald fc VennesI 1 1 99 ll ). It is important to recall that hydrogen can only accumulate in 
the photospheric regions, as a result of accretion, while heavier elements will slowly diffuse 
at the bottom of the helium conve ction zone on time scales of the order of 10^ years or less 
(see Table 1 of iDupuis et al.lll993l ). Since the bottom of the helium convection zone grows 
deeper with time (see Figure [3]), hydrogen becomes diluted in a larger volume. The measured 
hydrogen abundances in cool DB stars, and eventually DC white dwarfs, are thus governed 
by this balance between accretion and dilution within the helium convection zone. It is thus 
important to connect the constraints imposed by DBA stars with independent determina- 
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tions of the hydrogen abundances in cooler white dwarfs, such as DZA stars (jPufour et aL 



2007ah 



We show in Figure[57]the same results as in Figure [2^ but we added the H/He abundance 
ratios predicted from continuous accretion fro m the interstellar medium at various rates, from 



10 to 10 Mq yr ^ (see also Figure 10 of lVoss et al .112 00 71 for a sir nilar calculation). Als o 



shown are the hydrogen abundances for the DZA stars determined by lDufour et al.l (l2007al ). 
As can be seen, with the exception of the three hottest DBA stars, this range of accretion 
rates can easily account for the amount of hydrogen observed in DBA and DZA stars. 
However, some of the upper limits on the hydrogen abundance below ~ 20, 000 K are quite 
stringent, H/He < 10~^. These DB stars are observed in the same range as some DBA stars, 
characterized by hydrogen abundances of the order of H/He ~ 10~^, or a factor 100 higher. 
It is thus difficult to understand why accretion of hydrogen could be effective for some DBA 
stars, and not for other white dwarfs at the same temperature. The existence of such cool, 
hydrogen- deficient DB stars in our sample, and probably a significant fraction of cool DC 
and DZ stars as well, can only be interpreted as having evolved from hotter DB stars that 
contain negligible amounts of hydrogen in their atmospheres, without invoking any accretion 
mechanism whatsoever. 

If this interpretation is correct, the presence of hydrogen in DBA stars must be residual. 
Perhaps then, it is time to question the assumption of complete mixing of the hydrogen layer: 
in that case, it remains conceivable that the amount of hydrogen present is indeed of the 
order of ~ 10~^^ M©, i.e. the amount expected for the c onvective dilution of the hyd rogen 



atmosphere to occur below T^s ~ 20,000 K according to iMacDonald fc Venned (jl99ll ). but 
that it somehow fioats on top of the photosphere rather than being forcefully mixed by the 
helium convection zone. 

One particular object in our sample is LP 497-114 (1311-1-129), the DBA star with the 
largest hydrogen abundance, H/He ~ 10~^ (at Teg ~ 19, 000 K in Figure [25]). We have three 
independent spectroscopic observations for this object, one of which is a SDSS spectrum 
(which serves also as our Ha spectrum). Our best fit to these spectra are displayed in 
Figure [251 The two bottom fits are clearly at odds with what we can usually achieve in this 
temperature range, especially given the high signal-to-noise ratio of these observations. In 
particular, the region around He I A4471, as well as the core of most neutral helium lines, 
are poorly reproduced. If we consider only the blue spectra, we can achieve much better fits 
at Teff ~ 23, 000 K, but then the fine cores of Ha and He I A6678 are predicted much too 
shallow. Note that the fit to the SDSS spectrum (top) is less problematic, although the line 
cores are also predicted too shallow. Incidentally, LP 497-114 shows the worst agreement 
between optical and UV temperatures (Tuv ~ 27, 000 K in Figure [T5l). Perhaps the problem 
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with the modehng of this star is the assumption of a homogeneously mixed hydrogen and 
hehum atmospheric composition. This star might be in the process of being convectively 
mixed, with hydrogen constantly trying to float back to the surface. This could even explain 
the small spectroscopic variations from spectrum to spectrum observed in Figure If this 
is the case, the chemically homogeneous models used here no longer apply. This might also 
be the case, but to a lesser extent, with other DBA stars in our sample. 



5.5. The Instability Strip of the Pulsating V777 Her Stars 



Our first assessment of the instabi lity strip of the pu l sating DB (V777 Her) stars was 
presented in the spectroscopic study of iBeauchamp et al.l (119991 ). The picture at that time 
was complicated by the fact that only blue spectroscopic observations were available, and 
the insufficient constraints on t he hydrogen content is these stars did not allow us to choose 
the most appropriate solution. iHunter et al.l (l200l[ ) presented an update of the instability 
strip of the V777 Her stars using new spectroscopic observations at Ha, with the specific 
aim of constraining the hydrogen abundance in these stars. The results, shown in Figure 
2 of Hunter et al., reveal an instability strip that contains at least one non- variable white 
dwarf, as well as a significant number of unknowns. Since this last report, one obje ct in this 
sampl e, PG 2246+121, has been identified as a new variable DB white dwarf by [Handler 
J200l[ ). 



We show in Figure [29] the results from our revised spectroscopic analysis. The V777 Her 
stars are identified in Table 1. We distinguish in this plot variable and non- variable white 
dwarfs, but also DB and DBA stars. Above T^g ~ 20, 000 K, there are only 3 DBA stars 
in our sample with hydrogen abundances around H/He ~ 10~^'^, i.e. about one dex larger 
than the upper limits derived for DB stars in the same range of temperature, and all three 
are pulsating white dwarfs. Given our results presented so far, DB and DBA stars in this 
range of effective temperature seem to form two distinct classes of objects, with probably 
different progenitors: hot DB stars that once resided in the gap in one case, DA stars with 
thin hydrogen layers (Mh ~ 10~^^ Mq) in the other case. Our results displayed in Figure 
29] are consistent with this idea, in the sense that we can easily identify two distinct pure 
instability strips, within the uncertainties, one for the pure helium-atmosphere DB stars, 
and one for the DBA stars. 



Also displayed in Figure [29^ are the theoretica l blue edges for pure helium envelope mod- 
els repr oduced from Figure 3 of IBeauchamp et al.l (119991) . based on th e nonadiabatic calcula- 
tions of iFontaine fc Brassard] (11997] ) and lBrassard fc Fontaine] (119971 ). for various versions of 
the mixing-length theory. For DB stars, the observed blue edge near Teg ~ 30, 000 K suggests 
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a convective efficiency somewhere between the ML2 [a = 1.0) and MLS (= ML2/a = 2.0) 
versions of the mixing-length theory, a parameterization entirely consistent with that adopted 
in our model atmosphere calculations (ML2/a; = 1.25). The location of the red edge of 
the DB instability strip is located near T^s ~ 25, 000 K, for a corresponding strip about 
5000 K wide. Two objects, shown as black open circles, are located at the red edge — KUV 
03493+0131 (0349+015) and P G 2354+159 — but only t he latter has been confirmed as a 



non- variable, to our knowledge (IRobinson fc Wingetlll983[ ). 



The boundaries of the instability strip for the DBA stars are not as well defined as for 
the DB white dwarfs since we lack non-variable DBA stars both hotter and cooler than the 
DBA pulsators studied here, which would allow us to define empirically the edges of the 
instability strip (for instance, there is no DBA star in our sample between Tgfj ~ 23, 000 K 
and 20,000 K). Nevertheless, the existence of DBA pulsators cooler than the coolest variable 
DB stars suggests that the presence of small traces of hydrogen (H/He ~ 10~^) in the partial 
ionization zone would lower the temperature at which helium-atmosphere white dwarfs can 
pulsate. Interestingly, the effect produced by the presence of small amounts of hydrogen 
on the theoretical b l ue ed ge of the instability strip has been presented in Figure 10 of 



Fontaine fc BrassardI (120081 ). Indeed, the blue edge becomes cooler, but the effect is rather 
small (~ 500 K), although the authors considered a trace of only H/He = 10~^ in their 
calculations, and it is expected that the effect would be considerably more important with 
the higher hydrogen abundances inferred in our analysis. 

What finally comes out of these results is that DB white dwarfs cannot be considered 
as a homogeneous class of objects when studying the location of the instability strip because 
of this additional parameter, the hydrogen abundance. 



5.6. The Case of HE 2149-0516 



HE 2149—0516 is a DAB star discussed briefiy by IVoss et al.l (120071 ). who describe its 
spectrum as exhibiting strong Balmer lines and weaker but strong neutral helium lines. The 
authors had difficulty fitting this object with either pure hydrogen or pure helium models, 
and only a pure hydrogen fit to Ha suggested a temperature near Tes ~ 30, 000 K, with an 
extremely low surface gravity of log ~ 7. Our normalized spectrum for this peculiar white 
dwarf is displayed in Figure [201 The absolute energy distribution (not shown here) is almost 
fiat in the wavelength range displayed here, so this cannot be a hot white dwarf. Our best 
solution with our grid of DB models indeed yields T^s = 11,370 K, logg = 7.6, and log 
H/He = —4.5, but our spectroscopic fit completely fails to reproduce both the helium and 
hydrogen lines adequately. Instead, we could achieve a much better fit by assuming that HE 
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2149—0516 is an unresolved DA + DB (or perhaps DBA) composite system. In this case, each 
component of the system dilutes the absorption features of the other component. Given that 
there are simply too many free parameters to fit and that the observed features are relatively 
weak, we kept both values of log fixed in our fitting procedure, as well as the value of H/He 
for the DB star, and experimented with various values of these parameters. Our best fit is 
displayed in Figure [301 This is admittedly not a perfect fit, especially in the blue region of 
the spectrum, but given that both the DA and DB components fall in a temperature range 
which is problematic in terms of the modeling (i.e., the high- log (7 problem), we are rather 
satisfied with the quality of our fit. Note that the combined model fluxes also reproduce the 
slope of the observed energy distribution perfectly (not shown here). 



6. Conclusion 

We presented a comprehensive analysis of 108 DB white dwarfs using model atmosphere 
fits to blue and red spectroscopic observations. As in previous investigations, we showed that 
high signal-to-noise spectroscopic data near Ha are crucial to determine, or constrain, the 
hydrogen abundance in these white dwarfs. For the hottest stars in particular, such data 
allowed us to pinpoint more accurately the effective temperatures and surface gravities. 
While we demonstrated that these parameters depend on the assumed convective efficiency 
used in our model atmospheres, we also showed that a detailed parameterization of the 
mixing-length theory is more difficult to achieve than in DA stars, partly because of the lack 
of sensitivity of the UV energy distributions to this parameter. 

The mean mass of DB white dwarfs appears to be significantly larger than that of DA 
stars, especially at the cool end of the DB sequence. Since this corresponds to a temperature 
regime where the physics of line broadening becomes more questionable, van der Waals 
broadening in particular, we used trigonometric parallax measurements to demonstrate that 
the spectroscopic distances were in excellent agreement with those obtained from parallaxes, 
giving us confidence in our ability to model the DB stars with sufficient accuracy. We found 
no significant differences between the mass of DB and DBA stars, however, suggesting that 
both populations may have a common origin. DBA stars, on the other hand, are usually 
found at much lower temperatures than DB stars. While the luminosity functions of DA 
and DB stars identified in the PG survey indicate that 20% of all white dwarfs are DB stars, 
this appears to be true only below Tefj ~ 17, 000 K (i.e. Mboi > 9.5 in Figure [21]) . At higher 
temperatures, only ~ 9% of all white dwarfs are DB stars. This clearly indicates that the 
sudden increase in the ratio of DB to DA white dwarfs is the result of the transformation of 
DA stars into DB stars, but at a much lower temperature than the canonical value defined by 
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th e red edge of the DB ga p — or DB deficiency — near 30,000 K, a conclusion also reached 



by lEisenstein et al.l ( 120061 ) 



The global picture that ultimately emerges from our study is that there are at least 
two evolutionary channels feeding the DB white dwarf population. A small fraction of white 
dwarf stars are born as hydrogen-deficient stars, with negligible amounts of hydrogen in 
their envelope. These will remain hydrogen-deficient throughout their subsequent evolution, 
and will never become DA stars. The existence of extremely hot DB stars in the DB gap 
certainly supports this conclusion. We suggest that the hottest DB stars in our sample with 
Tefj > 20, 000 K descend directly from this evolutionary channel. In the bin centered on 
-/^boi = 7.0 in Figure 1211 we actually find only 3 DB stars, which represent 5.6% of the entire 
white dwarf population in this particular bin. If we apply the same fraction to the brighter 
bin centered on Mboi = 6.0 [T^s ~ 35,000 K), which contains 24 DA stars, we should only 
expect a single DB white dwarf in this particular bin. It is then perhaps not surprising that 
none were found in the PG survey. Note that a hot DB star in the gap was actually identified 
in the KUV survey, but this star was hidden in an unresolved DA + DB double degenerate 



system (ILimoges et al.ll2009l ). 



The hydrogen-deficient evolutionary channel is also supported by observations at low 
effective temperatures {T^s ~ 15, 000 K) where we see a significant number of DB stars with 
extremely low hydrogen abundances (H/He < 10~^). The mere existence of these almost 
pure helium atmospheres rules out the accretion from the interstellar medium scenario as 
the most likely source of hydrogen in the DBA stars found in the same temperature range. 
The only alternative, which is the second evolutionary channel, is that hydrogen must have 
a primordial origin, resulting from the transformation of a thin hydrogen-atmosphere DA 
white dwarf into a DB degenerate. The observed hydrogen abundances in DBA stars are 
too large, however, to be explained by a simple model where hydrogen is thoroughly mixed 
within the helium convection zone. Instead, we proposed a model where hydrogen tends to 
float on top of the photosphere rather than being completely mixed in the helium convective 
envelope. Since DBA stars appear in large numbers only below Teg ~ 20, 000 K, this suggests 
hydrogen layer masses of Mh ~ 10~^^ Mq for most DA progenitors, although the existence 
of a few hotter DBA stars may indicate even thinner hydrogen layers in some cases. 

In principle. Hot DQ stars must also be affecting the luminosity function of DB stars, 
either when DB stars are transformed into Hot DQ stars near Tcs ~ 24,000 K (i.e., the 
hottest DQ stars currently known), most likely through convective mixing, or when Hot DQ 
stars apparently return to being DB stars near T^s ~ 18, 000 K (i.e. the coolest DQ stars 
currently known), through a process currently unknown. But the number of known Hot DQ 
stars in the SDSS sample is so small, only 14 known to date, that the contribution of these 
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stars to the luminosity function of DB stars is most certainly negligible. 



We may speculate as to the origin of a larger mean mass for DB stars compared to DA 
stars. First, this may be an artifact associated with the model atmospheres. Indeed, the 
increase in mass becomes particularly important when the atmospheres of DB stars become 
strongly convective (see Figures [3] and [2T!) . Perhaps this increase in mass is analog to that 
observed in cool (Teg < 12, 000 K) DA stars, where the most likely explanation for this 
phenomenon has been attribute d to the limitations of the mixing-length th eory used in the 
model atmosphere calculations (IKoester et al.l l2009l : iTremblay et al.ll2010l ). It is possible 
that such limitations also apply to DB atmospheres. If the mean mass of DB stars is indeed 
significantly larger than that of DA stars, this implies that thei r progenitors on the mai n 
sequence must have been more massive as well. According to IWerner fc Herwid (120061 ) . 
the mechanism by which hydrogen is destroyed — during the post-asymptotic giant branch 
(AGB) phase — to the level required to account for the almost pure helium-atmosphere DB 
stars discussed here, or even the DBA stars in our sample, is probably caused by a very 
late helium-shell flash, or an AGB final thermal pulse. It is not inconceivable that such a 
mechanism is more efficient in more massive stars on the main sequence. 

Our next task will be to study the numerous DB white dwarfs discovered in the Sloan 
Digital Sky Survey, even though the average signal-to-noise ratio of the spectroscopic ob- 
servations is admittedly much lower than that of our sample (see Figure [6]). Future work 
should also include a proper treatment of evolutionary models with mixed hydrogen and 
helium abundances to evaluate better the total hydrogen content of DBA stars. 
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Table 1. Atmospheric Parameters of DB and DBA Stars 
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-4.84 


(1.41) 


0.63 


(0.03) 


10.57 


-1.43 


14.72 


67 


4.48(- 


-7) 


7.67 




1542-275 


LP 916-27 


12,700 


(385) 


9.13 


(0.15) 




-4.95 


(0.58) 


1.24 


(0.06) 


13.72 


-3.23 


15.49 


22 






9.19 




1545+244 


Ton 249 


12,840 


(332) 


8.18 


(0.13) 




-5.02 


(0.11) 


0.70 


(0.08) 


11.93 


-2.51 


15.78 


58 


2.88(- 


-6) 


8.63 




1551+175 


KUV 15519+1730 


15,550 


(378) 


7.96 


(0.11) 




-4.38 


(0.07) 


0.57 


(0.06) 


11.13 


-2.04 


17.50 


188 






8.24 


5 


1557+192 


KUV 15571+1913 


19,570 


(558) 


8.15 


(0.07) 




-4.37 


(0.30) 


0.68 


(0.04) 


10.92 


-1.75 


15.40 


78 






8.04 




1610+239 


PG 1610+239 


13,360 


(334) 


8.16 


(0.12) 


< 


-5.58 


(0.30) 


0.69 


(0.07) 


11.81 


-2.43 


15.34 


50 


2.42(- 


-6) 


8.57 




1612-111 


GD 198 


23,420 


(1781) 


7.96 


(0.06) 


< 


-4.75 


(2.11) 


0.59 


(0.04) 


10.44 


-1.32 


15.53 


104 






7.53 




1644+198 


PG 1644+199 


15,190 


(383) 


8.12 


(0.10) 


< 


-5.42 


(0.31) 


0.66 


(0.06) 


11.43 


-2.18 


15.20 


56 


1.41(- 


-6) 


8.38 




1645+325 


GD 358 


24,940 


(1115) 


7.92 


(0.06) 


< 


-4.58 


(0.88) 


0.57 


(0.03) 


10.33 


-1.19 


13.65 


46 


3.31(- 


-7) 


7.36 


3 


1654+160 


PG 1654+160 


29,410 


(1613) 


7.97 


(0.08) 


< 


-3.98 


(1.19) 


0.60 


(0.04) 


10.18 


-0.92 


16.55 


187 


2.62(- 


-7) 


7.07 


3 


1703+319 


PG 1703+319 


14,430 


(362) 


8.45 


(0.11) 




-5.51 


(0.34) 


0.88 


(0.07) 


12.09 


-2.48 


16.25 


67 


3.38(- 


-6) 


8.67 




1708-871 


L7-44 


23,980 


(1686) 


8.05 


(0.06) 


< 


-4.69 


(1.93) 


0.64 


(0.03) 


10.54 


-1.33 


14.38 


58 






7.55 




1709+230 


GD 205 


19,610 


(507) 


8.09 


(0.06) 




-4.00 


(0.13) 


0.65 


(0.04) 


10.81 


-1.71 


14.90 


65 






7.99 




1726-578 


L204-118 


14,320 


(341) 


8.20 


(0.08) 




-5.46 


(0.19) 


0.71 


(0.05) 


11.70 


-2.. 33 


15.27 


51 






8.51 




1822+410 


GD 378 


16,230 


(385) 


8.01 


(0.08) 




-4.45 


(0.06) 


0.60 


(0.05) 


11.09 


-2.00 


14.39 


45 






8.22 





Table 1 — Continued 



vv u 


JN cinic 


J-cB 




lo 


g 9 




log H/He 


M/Mq 


Mr . 


log 1^1 l^Q 


T/ 
V 




1/ Vmax 


log T 


NT /-if no 


1940+374 


L1573-31 


16,630 


(431) 


8.07 


(0.08) 


< 


-5.17 


(0.28) 


0.64 


(0.05) 


11.14 


-1.99 


14.51 


47 






8.23 




2034-532 


L279-25 


17,160 


(406) 


8.48 


(0.07) 


< 


-5.74 


(0.57) 


0.90 


(0.05) 


11.73 


-2.19 


14.46 


35 






8.49 




2058+342 


GD 392A 


12,220 


(423) 


9.09 


(0.20) 


< 


-5.54 


(2.56) 


1.23 


(0.08) 


13.72 


-3.26 


15.68 


24 






9.22 




2129+000 


G26-10 


14,380 


(351) 


8.26 


(0.14) 


< 


-6.48 


(1.66) 


0.75 


(0.09) 


11.79 


-2.36 


15.27 


49 






8.55 


2 


2130-047 


GD 233 


18,110 


(427) 


8.11 


(0.08) 




-5.76 


(1.32) 


0.66 


(0.05) 


11.01 


-1.86 


14.52 


50 






8.13 




2144-079 


G26-31 


16,340 


(410) 


8.18 


(0.07) 


< 


-6.22 


(1.45) 


0.70 


(0.04) 


11.36 


-2.09 


14.82 


49 






8.33 


2 


2147+280 


G188-27 


12,940 


(400) 


8.85 


(0.20) 


< 


-5.58 


(0.72) 


1.12 


(0.10) 


13.09 


-2.96 


14.68 


20 






9.12 




2222+683 


G241-6 


15,230 


(380) 


8.20 


(0.11) 


< 


-5.42 


(0.43) 


0.71 


(0.07) 


11.56 


-2.23 


15.65 


65 






8.43 


1 


2229+139 


PG 2229+139 


14,940 


(357) 


8.18 


(0.09) 




-4.73 


(0.07) 


0.70 


(0.05) 


11.57 


-2.25 


15.99 


76 


1.70(- 


6) 


8.44 




2234+064 


PG 2234+064 


23,770 


(1770) 


8.07 


(0.06) 


< 


-4.72 


(2.10) 


0.65 


(0.04) 


10.58 


-1.36 


16.03 


123 


4.51(- 


7) 


7.59 




2236+541 


KPD 2236+5410 


15,470 


(386) 


8.30 


(0.09) 


< 


-5.38 


(0.27) 


0.78 


(0.06) 


11.68 


-2.26 


16.19 


79 






8.48 




2246+120 


PG 2246+121 


27,070 


(1521) 


7.92 


(0.07) 


< 


-4.27 


(0.67) 


0.57 


(0.04) 


10.23 


-1.04 


16.73 


199 


2.86(- 


7) 


7.20 


3 


2250+746 


GD 554 


16,370 


(386) 


8.16 


(0.06) 


< 


-5.22 


(0.12) 


0.69 


(0.04) 


11.32 


-2.07 


16.69 


118 






8.31 




2253-062 


GD 243 


17,190 


(437) 


8.07 


(0.10) 




-4.33 


(0.12) 


0.64 


(0.06) 


11.05 


-1.93 


15.06 


63 






8.18 




2310+175 


KUV 23103+1736 


15,170 


(373) 


8.37 


(0.09) 


< 


-5.43 


(0.26) 


0.82 


(0.06) 


11.84 


-2.34 


15.88 


64 


2.40(- 


6) 


8.56 




2316-173 


G273-13 


12,610 


(424) 


9.11 


(0.19) 




-4.85 


(0.56) 


1.23 


(0.07) 


13.70 


-3.23 


14.08 


11 






9.19 




2328+510 


GD 406 


14,460 


(362) 


8.06 


(0.12) 


< 


-5.51 


(0.35) 


0.63 


(0.07) 


11.47 


-2.23 


15.09 


53 






8.42 




2354+159 


PG 2354+159 


24,830 


(1671) 


8.15 


(0.06) 


< 


-4.59 


(1.78) 


0.70 


(0.04) 


10.67 


-1.34 


15.78 


105 


5.00(- 


7) 


7.57 





Note. — (1) Solution with log Ca/He = —7.0; (2) Limits on the hydrogen abundance based on the absence of Hci from Table 1 of lVoss et aI.l[2007l L (3) Variable 
white dwarf of the V777 Her class; (4) PG star not in the complete sample; (5) Solution with log Ca/He = —6.5. 
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Table 2. Parameters of DB Stars with Trigonometric Parallaxes 







7r 


0"7r 




^spcc 










WD 


Name 


(mas) 


(mas) 


(pc) 


(pc) 


(K) 


log g 


M/Mq 


Notes 


0002+729 


GD 408 


28.80 


(4.7) 


34.7 


32.2 


14,410 


8.26 


0.75 


1 


0017+136 


Feige 4 


30.00 


(9.4) 


33.3 


75.9 


18,130 


8.09 


0.65 


1 


0615-591 


L182-61 


27.50 


(1.0) 


36.4 


37.5 


15,750 


8.04 


0.61 


2 


1333+487 


GD 325 


28.60 


(3.2) 


35.0 


35.3 


15,320 


8.03 


0.61 


1 


1425+540 


G200-39 


17.30 


(3.9) 


57.8 


56.2 


14,490 


7.95 


0.56 


1 


1645+325 


GD 358 


27.30 


(3.3) 


36.6 


46.0 


24,940 


7.92 


0.57 


1 


1940+374 


L1573-31 


20.30 


(2.8) 


49.3 


47.2 


16,630 


8.07 


0.64 


1 


2129+000 


G26-10 


20.26 


(2.0) 


49.4 


49.6 


14,380 


8.26 


0.75 


2 


2144-079 


G26-31 


14.40 


(5.8) 


69.4 


49.2 


16,340 


8.18 


0.70 


1 


2147+280 


G188-27 


28.30 


(3.0) 


35.3 


20.8 


12,940 


8.85 


1.12 


1 


2222+683 


G241-6 


14.70 


(4.4) 


68.0 


65.8 


15,230 


8.20 


0.71 


3 



Note. — (1) YPC, (2) iGould fc Chanamell2004l. C3'l lDahn et"al]|l988l . 
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Table 3. White Dwarf Space Density^ in the PG Survey 



Mboi DB 




DA 




Total 




6 — 




2.43 X 10" 


-6 


2.43 X 10" 


-6 


7 7.28 X 10- 


-7 


1.23 X 10- 


-5 


1.30 X 10- 


-5 


8 3.10 X 10- 


-6 


2.48 X 10- 


-5 


2.79 X 10- 


-5 


9 4.57 X 10- 


-6 


5.15 X 10" 


-5 


5.61 X 10- 


-5 


10 2.57 X lo- 


-5 


6.00 X 10- 


-5 


8.58 X 10- 


-5 


ll 1.73 X 10- 


-5 


6.11 X 10- 


-5 


7.85 X 10- 


-5 


^AU units in pc" 


-3 


mag-^. 
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Fig. 1. — Left panel: synthetic spectra of pure helium DB models [logg = 8 and ML2/a = 
1.25) at various effective temperatures (blue) normalized to a continuum set to unity and 
offset from each other by a factor of 0.4 for clarity; each model spectrum is compared to 
a spectrum at T^s = 24,000 K (red). Middle panel: same as left panel but for synthetic 
spectra at \ogg = 7.5 (red) and logg = 8.5 (blue). Right panel: same as left panels but for 
synthetic spectra calculated with ML2/q; = 0.75 (red) and a = 1.75 (blue). 

Fig. 2. — Variation of the equivalent width of He I A4471 (measured from 4220 to 4625 A) 
as a function of effective temperature for our grid of model spectra with various hydrogen 
abundances (left panel), surface gravities (middle panel), and convective efficiencies (right 
panel). Also shown as red circles are the measured equivalent widths in our sample of DB 
stars, arbitrarily located at Tes = 22, 000 K. 

Fig. 3. — Location of the helium convection zone (hatched region) as a function of effective 
temperature in the pure helium envelope of a 0.6 Mq DB white dwarf calculated with the 
ML2/a = 0.6 version of the mixing-length theory (from G. Fontaine & P. Brassard 2006, 
private communication). The depth is expressed as the fractional mass above the point of 
interest with respect to the total mass of the star. The thick solid line corresponds to the 
photosphere (r/j ~ 1). 

Fig. 4. — Distribution of DB stars in our sample as a function of the V magnitude. 

Fig. 5. — Optical (blue) spectra for our complete sample of DB stars. The spectra are 
normalized at 4500 A and shifted vertically from each other by a factor of 0.5 for clarity. 
The effective temperature decreases from upper left to bottom right. 

Fig. 6. — Distribution of signal-to-noise ratios for the 108 optical spectra secured for the 
analysis of the DB stars in our sample (hatched histogram). Also shown are the corre- 
sponding values for 744 DB stars identified in the Data Release 4 of the SDSS (thick solid 
line) . 

Fig. 7. — Misclassified spectra in our original target list of DB white dwarfs. 

Fig. 8. — Optical spectra in the red for the DB stars in our sample, when available, ordered 
in right ascension. The spectra are normalized at 6200 A and shifted vertically from each 
other by a factor of 0.5 for clarity. The dotted line indicates the location of Ha. Spectra 
labeled in red are taken from the SDS S. Ha spect r a from the SPY survey, also used in our 



analysis of DBA stars, are displayed in IVoss et al.l (120071 ) 



Fig. 9. — Example of the procedure used to define the continuum. In the top panel, the 
observed spectrum is fitted with a model spectrum multiplied by a high order polynomial (up 
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to A^) to achieve the best possible match (shown in red); in this case the resulting atmospheric 
parameters are meaningless since too many free parameters are used in the fitting procedure. 
The continuum flux is then defined by this fitting function at some predefined wavelength 
points (shown as blue tick marks in the bottom panel), and these are then used to normalize 
each segment of the spectrum to a continuum set to unity, as shown in the bottom panel. 

Fig. 10. — Examples of spectroscopic fits for two DB white dwarfs using a cool (bottom 
panels) and a hot (top panels) initial estimate of the effective temperature in our fitting 
procedure. 

Fig. 11. — Example of a full spectroscopic fit where the Ha line profile, shown in the insert, 
is used to measure, or constrain, the hydrogen abundance of the overall solution. 

Fig. 12. — Limits on the hydrogen abundance set by our spectroscopic observations at Ha 
and H/3. We estimate the limit of detectability at an equivalent width of 200 mA and 300 
mA, respectively. 

Fig. 13. — Our best fits to the blue spectra of the strongest DBZ stars in our sample. The 
region near Ha is not displayed here. 

Fig. 14. — Distribution of surface gravity as a function of effective temperature for all DB 
white dwarfs in our sample for various parameterization of the convective efficiency. Also 
shown in each panel is a 0.6 Mq evolutionary sequence. 

Fig. 15. — Top two panels: Comparison of optical and UV temperatures for 34 DB stars in 
our sample with available /[/i? spectra, using two differen t parameterization of the convective 



efficiency. Bottom panel: UV temperatures obtained by ICastanheira et al.l (l2006[ l based on 



ML2/a = 0.6 models compared with our own temperature estimates based on similar models. 

Fig. 16. — Comparison of UV temperatures obtained from model spectra calculated with 
different convective efficiencies (top panel), different hydrogen abundances (middle panel), 
and different surface gravities (bottom panel). In each panel, the UV temperatures on the 
y-axis, obtained from models with parameters given in the legend, are compared against UV 
temperatures obtained from pure helium model atmospheres at log (7 = 8 and ML2/a = 1.25. 

Fig. 17. — Comparison of Tefj and logg determinations for 28 DB stars in our sample with 
mutliple spectroscopic observations. 

Fig. 18. — Sample fits for several DB and DBA stars in our sample; the atmospheric pa- 
rameters (Teff, logg, and log H/He) of each object are given in the figure. The region near 
Ha (right panel) is used to measure, or constrain, the hydrogen abundance. In the case of 
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DB stars, these high signal-to-noise spectra provide upper limits on the hydrogen-to-hehum 
abundance ratio. 

Fig. 19-^ Comp arison of T^s and logg measurements between our analysis and that of 
Voss et al.l ( 120071 ) as a function of our T^s values (this work — Voss in each case). The 
dashed lines represent a perfect agreement. Objects shown as open circles correspond to 
cool DB white dwarfs for which Voss et al. assumed a value of logg = 8. The two objects 
labeled in the figure correspond to (1) PG 1115+128 and (2) KUV 03493+0131 (0349+015). 



Fig. 20. — Spectroscopic solutions obtained by I Voss et al.l ( 120071 ) for PG 1115+158 and KUV 



03493+0131. These should be contrasted with our solutions displayed in Figures [TTI and [TS| 
respectively. 

Fig. 21. — Distribution of mass as a function of effective temperature for the 61 DB (red 
symbols) and 47 DBA (blue symbols) stars in our sample. Also shown as dotted lines are the 
theoretical isochrones from our evolutionary models, labeled as logr where r is the white 
dwarf cooling age in years. 

Fig. 22. — Comparison of distances inferred from trigonometric parallax measurements (-Dtt) 
with those derived from spectroscopy (Dgpec) taken from Table 2. The objects labeled in 
the figure correspond to (1) Feige 4 (0017+136), (2) CD 358 (1645+325), and (3) G188-27 
(2147+280). 

Fig. 23. — Mass distribution of all DB and DBA white dwarfs in our sample (hatched 
histogram). If the most massive DB stars in this distribution are excluded (see text), the 
mean mass becomes of (M) = 0.671 Mq with a standard deviation of um = 0.091 Mq. Also 
shown as a thick solid line is the mass distribution of the DA stars in the PG sample. 

Fig. 24. — Luminosity function derived for all DB and DA stars from the complete PG 
sample presente d in bolometric mag nitude bins, assuming a scale height for the Galaxy of 



zo = 250 pc (see lLiebert et al.l l2005l for details). 



Fig. 25. — Hydrogen-to-helium abundance ratio as a function of effective temperature for 
all DB (white symbols) and DBA (red symbols) white dwarfs in our sample. Limits on the 
hydrogen abundance set by our spectroscopic observations at Ha (lower dotted line) and H/3 
(upper dotted line) are reproduced from Figure |T2j The hydrogen abundances for DB stars 
represent only upper limits. 

Fig. 26. — Total hydrogen mass as a function of effective temperature for all DB (white 
symbols) and DBA (red symbols) in our sample. The values for DB stars represent only 
upper limits. 
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Fig. 27 . — Same as Figure |25| but including also the results for DZA stars from lDufour et al. 



( l2007a[ ). Solid lines represent the expected abundances for continuous accretion of material 



from the interstellar medium with accretion rates of 10 to 10 ^'^ Mq yr ^. 

Fig. 28. — Our best fit to three independent spectroscopic observations of LP 497-114 
(1311+129), the DBA star with the largest hydrogen abundance in our sample. Our adopted 
solution is derived from the bottom spectrum, obtained in 2010. The other spectra have been 
secured in 1995 (middle) and 2004 (top; from SDSS). 

Fig. 29. — Location of the V777 Her stars in a mass versus effective temperature diagram. 
Variable DB white dwarfs are represented by filled symbols while open symbols correspond 
to photometrically constant, or unknown, objects. DB and DBA stars are shown in black 
and red, respectively. The dashed lines correspond to theoretical blue edges for pure helium 
envelope models with convective efficiencies given by, from left to right, ML3, ML2, and 
MLl. 

Fig. 30. — Our best fit to HE 2149—0516 assuming composite DA + DB models. The 
atmospheric parameters of each component are given in the figure; the log (7 values have 
been kept constant in the fitting procedure, and a hydrogen abundance of H/He = 10~^ was 
also assumed for the DB component. Note that the relative fiux scale starts at 0.4 so the 
absorption lines are particularly weak as a result of the dilution by each component of the 
system. 
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Figure [3] 
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Figure [6] 



- 51 - 






0447+176 



1053-092 



1302+317 



1614+270 



1632+222 




1816+667 



2006+481 




2146+412 




2234+033 




2358+155 



_±_ 



4000 



4500 



5000 



Wavelength (A) 



Figure [7] 



- 52 - 




o in 



Figure [8] 




Wavelength (A) 




Figure [10] 



- 55 - 



PG 1115+158 




log H/He = -3.84 F i " 

6400 6600 6800 



3800 4000 4200 4400 4600 4800 5000 

Wavelength (A) 



Figure [TT] 
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